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Integral sliding mode controller NPC is known by its fast dynamic and high capacity to reject disturbances. 
MPPT perturb and Observe The hybrid system is modelled in MATLAB/Simulink software. During 
Non-linear predictive control simulations, the DISMC-MPPT is compared with other techniques such as 


sliding mode controller (SMC) MPPT and integral SMC MPPT, the DISMC 
provides the best tracking performances under different irradiances. 
Moreover, the designed controller for the bidirectional converter regulates 
the DC-link voltage with better performances than the classical PI controller. 
Lastly, the NPC regulates the speed of the IM with high robustness. 
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1. INTRODUCTION 

The sliding mode control (SMC) is suitable for fast maximum power point tracking (MPPT) because 
it is known by stability, fast-response and non-sensitivity to parameter variation [1]. The conventional 
hysteresis-modulation (HM) based SMC suffers from variable switching frequency and from the chattering 
phenomenon due to the high switching frequency operation. These drawbacks motivated the researchers to 
design an indirect SMC based on the pulse width modulation (PWM) technique. The conventional SMC with 
PWM exhibits an unwanted steady-state error and slow response [2], to improve the response of indirect 
SMC an integral term is added to the existing sliding surface to constitute the integral sliding mode controller 
(ISMC) [3, 4]. However, the construction of the indirect form of the ISMC derivates the state variables of the 
switching surface (S=0), thus, the variable f x; dt disappears from the equivalent control term ueq, and the 
correction of steady-state error is deteriorated [5, 6]. Therefore, an integral term is added for the second time 
J (Cf x; dt))dt to nullify the steady-state error introduced by the indirect ISMC, this new technique is called 
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PWM based double integral SMC (DISMC) [2]. The DISMC shows a fast response and reduces the steady- 
state error [4]. The bidirectional (BDC) DC-DC converter is the key factor of energy management between 
PV source, battery, and different loads [7, 8]. However, the use of two PI controllers to control the 
bidirectional converter does not guarantee good performances when the system operates outside of the 
operating point [9, 10]. For that reason, the features of the ISMC theory are exploited to control the 
bidirectional converter in this paper. 

The discovering of field-oriented control (FOC) in 1972 has revolutionized the control theory of the 
induction motor (IM) [11-13]. However, the FOC technique uses coordinate transformation which makes it 
complex to implement, in addition, it is sensitive to internal parameters variation of the machine [14, 15]. To 
minimize the complexity of FOC and improving its dynamic new control techniques have been founded, 
namely sliding mode control (SMC), and direct torque control (DTC) [16, 17]. The IM model is non-linear in 
nature, furthermore, the IM is applied to external disturbances (load torque) and to internal parameter 
variations. Therefore, the robust non-linear predictive control (NPC) is applied to control the IM in this work. 
The NPC uses high order calculations to optimize a cost function; the complexity of calculation could be 
reduced using the Taylor series expansion of the output’s variables [18, 19]. 

This contribution proposes a high-performance control scheme for a hybrid water pumping system. 
Section 2 presents a brief description of the different stages of the pumping system. Section 3 presents the 
different control strategies for the pumping system. First, a cascaded controller based on DISMC is presented 
for MPPT. Afterward, two other controllers based on ISMC theory are designed to control the bidirectional 
flow of the energy. Finally, the control scheme for the IM based on NPC is explained. Simulation results and 
analysis are reported to section 4, followed by general conclusions. 


2. CIRCUIT CONFIGURATION 

A hybrid water pumping system is presented in Figure 1. The overall system comprises a PV array 
as a principal source of energy and a battery pack considered as a second power source. The system includes 
static DC-DC converters such as a unidirectional boost converter used for MPPT and a bidirectional buck- 
boost converter which ensures the bidirectional flow of the energy between the DC-bus and the battery. 
Finally, a two-level inverter controls an induction motor with a centrifugal pump. 
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Figure 1. Block diagram of the hybrid system 
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3. CONTROL STRATEGIES 
Diverse controllers are proposed to control the hybrid system, and are listed as follows: 
e A DISMC based on PWM is designed to track the maximum power point (MPP) with high efficiency. 
e The bidirectional buck-boost converter is controlled with the ISMC theory to regulate the DC-bus voltage. 
e Control of the induction motor with the generalized predictive controller. 


3.1. Maximum power point tracking control 

Figure 2 shows the MPPT control scheme, in which a perturb and observe (P&O) algorithm 
generates the reference voltage Vpy. Then, a DISMC is used to force the state trajectory of the boost 
converter to follow the reference voltage given by the P&O algorithm. 
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Figure 2. DISMC- MPPT algorithm control scheme 
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3.1.1. Design of a double integral sliding mode controller for MPPT 
Considering that the boost converter operates in continuous conduction mode (CCM), The average 
state-space model of the boost converter is presented in term of the switching signal u; as, 


x, = f(%)+ g(x )u; (1) 








L PG. 


py “in 


i T 
: T = Voy de K Ly — Vie 
Where, x, = i, Viv | , f(x) | Z 4 and g(x) = % | 
Where, rpy is the dynamic resistance of the PV array. The general control law to control the voltage vyv is: 


0 when S >0 
U, = 
'}1 when S <0 a 


Where S; is the switching surface, which is expressed as, 
S, =a, + a,€, Fae tae, (3) 
The terms a- a4 denote the sliding surface parameters, and e1- e4 are the error signals. 
e,=i,-i, / i,=A(V,-y,,) 
e, = Vi, = 
e = |V,- fv, at (4) 
e, = f| [(V;, - Av, at ar 
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Where, A is the amplified gain of the voltage error. The application of the indirect SMC is achieved through 
the computation of the equivalent control Ujeg. This can be derived from the invariance condition, S;=0, 





Si =a, €1+ a, €2+ a, e3t+ a, e4 =0 (5) 
ic, Vy e e 
Then, a O ey E (6) 
Vi Vi Yi Vie 
APL apL aL a,L 
With, A, =(——+ Z ), A= and A, =—— 
C, aC, a a 


Ai, Az, and A3 parameters are empirically determined according to the existence and stability criterions [2, 
4]. The control law of the indirect sliding mode control is derived using PWM technique, comparing control 
signal Vcontroi With a ramp signal Vrampı. 


pe = (VaT Vy) 7 Alc, T Ae, + A,e, (7) 


Vapi E p Yi 


3.2. The control strategy of the bidirectional converter 

The BDC is considered as the intermediate link between the battery and the DC bus. If the DC-link 
voltage is greater than the reference voltage Vac> Vac , the BDC operates in buck mode to store the excess of 
energy in the battery. On the other hand, if vac< Vac the battery discharges to inject the required current into 
the DC-bus, in this case, the BDC works in boost mode [20, 21]. The direction of the current is reversed from 
one mode to another. To control the BDC in both modes of operation two strategies of control are proposed 
as presented in Figure 3. 
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Figure 3. The proposed control scheme for the bidirectional converter 


3.2.1. Design of a double integral sliding mode controller for the boost converter mode 

Considering that the boost converter operates in CCM, and the chosen state variables are the 
inductor current itg and the DC-link voltage Vac. The dynamic of the boost converter is described with the 
differential equations in (8), where uz refers to the state of the switch Kg. 
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Xa = f(x,)+ g(x, )u, (8) 
; T , T 
— y i l y j; y : —] 
Where, X, = [in Val f(%)= a fate and &(X,) |i de 


Tac 18 the dynamic resistance of the DC bus. 
The sliding surface is formed with the state variables errors (9). 


S, =be, +b,e, +0,e, + b,e, (9) 
where e1-e4 are the current and voltage errors, and B is the gain which amplifies the voltage error. 


x 


e =i —lyp / is = BV; = Pvg) 
ers a — By, 
e; = | (e, +e, )at (10) 


e,= ffe +e,)dt |ar 


Substituting the dynamic of the step-up converter in the derivative of the sliding surface (S$2=0) gives, 


Vg UVa 





y IC IC 
S2 =b(# -BTH -b PT + by(e, +e) +b, |e +e dt (11) 


B dc dc 


The equivalent control signal U2eq of the boost converter is obtained by solving the following equation So=0, 


ly vy, bL, Di eae bL, 
u, =l- B+ Epy. (VV. - By +8 |V- By, jdt +711 BV, - By t 
eq2 Ci oe ite v, b, 3 de - pv,,) Vie ia 4! B dc — pva) dc -lig is |+ b, A K dc p a b, 7 [fB dc api) de — lip i d 


(12) 


The control law is derived using the PWM technique, comparing the voltage control signal Vcontroi2 With a 
ramp signal Vramp2. 


Vono = B (V; z Pid) a B, [BV z pv ,.) = bi + B[V; z pv, dt + B, [| BV; p Pv.) F lip Jar — Bic, T G, (Vs Vg) 
=G V; 


Vrapi 


(13) 


Here, the factor Gs=ß ( 0<G,<1) is used to downscale voltage magnitude to a practical level. B1-B3 
are constant parameters determined according to reachability and stability conditions [5, 22]. The duty ratio 
is multiplied with a signal of the pulse generator to ensure that the duty ratio is always less than 1. Figure 3 
illustrates the control scheme of the implemented DISMC-PWM. 
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3.2.2. Design of Integral sliding mode controller for the buck converter mode 

Considering that the buck converter operates in CCM, and the chosen state variables are the inductor 
current Ig and *° DC-link voltage Vac.The dynamic of the buck converter is described with the differential 
equations (14), where us refers to the state of the switch Kpp. 


x3 = f (x3) + g(a, )u, (14) 
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T T 
j 7 “vz, i v —1 
where, Xa =|lp V b : x,) =| —* de and x) =| LB 
s= lie Val . FC 3) z% i g(x; ) E | 
the switching surface S; is formed with the DC-link voltage error (e1), the integration of this error (e2) and the 


inductor current ip. The sliding surface S3 is expressed in (15), 
S = Ce +C, +i, (15) 


where, cı and cz are the sliding surface parameters determined according to reachability condition [9, 25], the 
errors are given by, 


x 
e =V Ya 
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+ 16 
e, = | Vi -v,,)at m 


the control discrete function of the buck converter is concluded from the transversality condition [25]. 
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the derivative of the switching surface is given in equation (18), 
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3.3. Non-linear predictive control of induction motor 
Figure 4 presents the non-linear predictive control of the induction motor. The mathematical model 
of the IM is presented in the two-dimensional stator reference frame (a-B) [26] is: 


x= f(x) +9(x)ult) (20) 
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where, Usa, Usp, Isa, isp are the stator voltages and the stator currents, respectively, and Pra, Orp are rotor fluxes. 
The constant parameters of the model are defined by, 








i L, 1 La 
o=1-| — |; K= f= hak) 
LL, oL,L, oL, L 


y 


where, Ls, L; are stator and rotor inductances, and Lm is the mutual inductance, Rs, R; are stator and rotor 
resistances, and T,= L,/R,; is rotor time constant, p number of poles pair, Tı load torque, J inertia coefficient, f; 
friction coefficient. The square of rotor flux and the speed are two chosen output to be regulated, 


y=h(x)=(y, y,) =(h(x) h) =(0, 67=6, thp ) (21) 


the controlled variables can be expressed with Lie derivatives as, 


YO! f hw 054 
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where, G; (x) = , L hEHL h OL, h oO, 0<i<2 


where, L; h j (x) operator is Lie derivative notation of the function h; with the respect to f(x). 
The basic idea behind the predictive control consists of the construction of a control law u(t) able to force the 


system trajectory to follow the desired trajectory in a future horizon (t+t;). The control law is obtained 
through the optimization of the cost function expressed as, 


J == [t+0-y, D ltt 2)-y, (tee) de (23) 


using Taylor series expansion of the outputs and of the reference outputs the cost function can be rewritten as 
[15, 18]: 


J = (YO -Y, © THY -Y, ©) (24) 
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the optimal control law is deduced by setting the following condition, 


0S GO [T T |M +a TG, ou = 0 (25) 
u 


finally, the control law u(t) is obtained by, 
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u(t) = -G (x) T, ie, Lo |m (26) 


the conditions { T +o + 0} and {Paro #() } should be fulfilled to inverse the matrix G;(x) [15]. 
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4. SIMULATION RESULTS AND ANALYSIS 

The hybrid pumping system has been tested in MATLAB / SimulinkTM software. The principal 
data of the hybrid pumping system are listed in Table 2. Figure 5 (a) shows that linear and sudden irradiances 
are applied to the system, and the temperature is fixed at 25 ° C value due to its little effect on power 
variation. Figure 5 (b) describes the evolution of PV voltage, the controller based on DISCM-MPPT tracks 
fastly the reference voltage Vpv provided by P&O algorithm. In addition, a zoomed view on the PV voltage 
curve illustrates that the voltage fluctuations are small when the MPP is reached. In order to make a fair 
comparison possible, the designed DISMC-MPPT is compared with both the conventional SMC-MPPT and 
the ISMC-MPPT under the same operating conditions of the system. Figure 5 (c) and Table 1 clearly show 
the improvement of the extracted power. From simulation data, it can be observed that the conventional SMC- 
MPPT presents the highest power fluctuations, the slower tracking speed and the largest steady-state error. 
The ISMC-MPPT increases the tracking speed, but the amplitude of power oscillation is not reduced 
significantly. In contrast, the DISM-MPPT reduces the amplitude of power oscillation and improves the other 
performances in comparison with ISMC-MPPT. 


Table 1. Comparison of the MPPT properties in terms of power oscillations and tracking speed during sudden 
change of irradiance 
Irradiance (W/m?) 


600 400 200 300 500 700 
MPPT Track. Power Track. Power Track. Power Track. Power Track. Power Track. Power 
Technique speed Oscill. speed Oscill. speed Oscill. speed Oscill. speed Oscill. speed Oscill. 
(ms) (W) (ms) (W) (ms) (W) (ms) (W) (ms) (W) (ms) (W) 
SMC 400 0.8 2 0.3 - 0.09 - 0.2 2 0.9 2 0.8 
ISMC 69 0.6 1.2 0.3 - 0.07 -0.16 1.85 0.9 1.8 0.8 
DISMC 69 0.5 1 0.2 - 0.05 - 0.1 1.5 0.3 1.5 0.6 


On the other side, the sliding mode control scheme for the bidirectional converter is compared with 
the conventional PI control scheme. Figures 6 (a) presents the regulated DC-link voltage with both control 
strategies. For the PI-based controller, it can be seen that the DC-link voltage moves away from its reference 
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when the bidirectional converter switches from one mode (charging/discharging) to another. On the other 
hand, the DC-link voltage is accurately regulated with the sliding mode controllers, in this case, the system 
presents an insignificant drop of voltage and short setting time when the operating mode changes. 
Furthermore, the voltage ripples are reduced with the proposed controller. The state of charge (SOC) shown 
in Figure 6 (b) reflects the battery’s charging/discharging modes. The SOC increases when the battery is 
charging and vice versa. Figure 6 (c) shows that the inductor Lg operates in CCM mode, the battery receives 
the current from the DC-bus during charging (ig negative) and delivers the current to DC-bus for the time of 
discharge (itp positive). 

It can be observed from Figure 7 (a) that the non-linear predictive controller tracks the reference 
speed with a very fast dynamic. The electromagnetic torque curve presents reduced ripples because the 
curves of stator currents are quite close to the sinusoid form, which is depicted in Figure 7 (b) and Figure 7 
(c). Figure 7 (d) shows that the decoupling between the flux and the torque is achieved accurately because the 
flux stays close to its reference. Since the NPC controller tracks the reference speed with high performance, 
the centrifugal pump parameters (flowrate, hydraulic power) are improved. Figure 7 (e) and Figure 7 (f) 
illustrate that these parameters follow exactly their references. 
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Figure 5. Simulation results of the PV array parameters under linear and sudden irradiances. 
(a) Irradiance, (b) PV voltage, (c) The extracted PV power for different methods 


Table 2. Hybrid water pumping system parameters 


Parameter Value 
Current at maximum power (Imp) 17.68A 
PV array Voltage at maximum power (V mp) 204V 
Maximum power (Pm) at STC 3606.72W 
Rated power (Pm) 2200W 
Rated speed (@,) 1430 rpm 
Induction Motor Drive Rated frequency (f) 50Hz 
Nominal voltage (L-L) 230V 
Rated flux (0,) 0.6Wb 
Rated pumping flow (Q,) 10m°*/h 
Cent tugat Pomp Rated pumping head (H,) 50m 
Pump constant (Kpump) 6.55x104W/(r/s} 
Battery Pack Nominal voltage (Vs) 192V 
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Figure 6. Simulation results of battery parameters under different load variations. 
(a) DC-link voltage for both PI and SMC methods, (b) state of charge, (c) the inductor current itg 
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Figure 7. Simulation results of the motor-pump parameters 
(a) rotor speed, (b) electromagnetic torque, (c) stator currents, (d) rotor flux, (e) pump flowrate, 
(f) hydraulic power 


5. CONCLUSION 

A robust control scheme of a hybrid water pumping system based on an induction motor was 
presented in this paper. First, a cascaded controller based on P&O and the indirect double integral SMC is 
designed to track the MPP. Seconde, the integral sliding mode control theory has been used to control the 
bidirectional buck-boost converter. The third stage of the system consists of a three-phase inverter which 
controls the speed and torque of an induction motor with the non-linear predictive control technique. 
Analyzing simulation results, it is found that the DISMC-MPPT tracks the MPP with better performances 
than both SMC-MPPT and ISMC-MPPT methods. The bidirectional converter controlled with the 
conventional PI controller is not able to preserve the desired response. However, the control scheme based on 
ISMC theory showed high robustness and operated properly in different irradiances and speeds. Moreover, 
the NPC exhibited high tracking performances of the speed, fast torque response and fewer ripples in the 
torque. Since the IM is controlled with high performances, the pump parameters (hydraulic power, flowrate) 
are improved. 
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